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In the study of weakly turbulent wave systems possessing incomplete self-similarity it is possible 
to use dimensional arguments to derive the scaling exponents of the Kolmogorov-Zakharov spectra, 
provided the order of the resonant wave interactions responsible for nonlinear energy transfer is 
known. Furthermore one can easily derive conditions for the breakdown of the weak turbulence 
approximation. It is found that for incompletely self-similar systems dominated by three wave 
interactions, the weak turbulence approximation usually cannot break down at small scales. It 
follows that such systems cannot exhibit small scale intermittency. For systems dominated by four 
wave interactions, the incomplete self-similarity property implies that the scaling of the interaction 
coefficient depends only on the physical dimension of the system. These results are used to build 
a complete picture of the scaling properties of the surface wave problem where both gravity and 
surface tension play a role. We argue that, for large values of the energy flux, there should be two 
weakly turbulent scaling regions matched together via a region of strongly nonlinear turbulence. 

PACS numbers: 04.30.Nk, 92.10.Hm, 92.10.Cg 



I. INTRODUCTION 



The time evolution of the average spectral wave-action density, rik, of an ensemble of weakly interacting dispersive 
waves is governed by the so-called kinetic equation. For a system with dispersion law Wk, dominated by 3-wave 
interactions with interaction coefficient ivkkik2i the kinetic equation is : 
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More generally, if the dominant nonlinear interaction is iV-wave, the schematic form of the kinetic equation is 

jL%n^-U{c.^)6{k){dkr-\ (3) 

Let us consider homogeneous, isotropic systems in physical dimension d. Let us further assume scale invariance with 
tJk and ijv homogeneous functions of their arguments. Denote the degrees of homogeneity of the dispersion, o^k, and 
the iV-wave interaction coefficient, Ljv, by a and respectively. Under these assumptions, the kinetic equation 
possesses exact stationary solutions, found originally by Zakharov in the early 70's, which carry constant fluxes of 
conserved quantities, such as energy or wave-action. These solutions are called Kolmogorov-Zakharov (KZ) spectra. 
The 3-wave kinetic equation admits a single KZ spectrum carrying a constant flux, P, of energy : 



nk 



= c(3)pU-('^-^+'". (4) 



The 4-wave kinetic equation conserves wave-action in addition to energy and thus admits a pair of KZ spectra, one 
carrying an energy flux, P, the other carrying a wave-action flux, Q. They are 

nk = cJ^VH-^(2^*+3'*), (5) 



The dimensional constants, C'^-', can be explicitly calculated. 
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Suppose we know that the physical system under consideration depends on only one dimensional constant. Such 
a wave system is said to possess incomplete self-similarity (ISS). Zakharov, Lvov and Falkovich have pointed out 0, 
chap. 3] that for such systems the scaling of the KZ spectra can be obtained from a dimensional argument. The 
dimensional argument for ISS systems uses only the scaling of the dispersion relation and not the scaling of the 
interaction coefRcients required for the more general scale invariant systems considered above. This fact has not been 
fully appreciated and is rarely used, despite the fact that most of the known wave turbulence systems are ISS, as will 
be shown in this paper. 

The dimensional analysis determines the scaling of the interaction coefficients for systems possessing ISS, a point 
that has been mostly overlooked before. The fact that the scaling exponent of the interaction coefficient is not 
an independent quantity may have some consequences for the practical applicability of some theoretical results on 
weak turbulence, where the scaling exponents of the interaction coefficients are regarded as arbitrary. In particular, 
Biven, Nazarenko and Newell0,|^ have recently pointed out that the weak turbulence approximation is almost never 
uniformly valid in k but rather breaks down either at large or small scales. The breakdown of weak turbulence at 
small scales is presumed to signal the onset of small scale intermittency. It is possible to use a simple dimensional 
argument to recover the criteria obtained in in the ISS case. One finds that the condition for breakdown at small 
scales is inconsistent for three wave systems. As a result one would not expect such systems to exhibit small scale 
intermittency. 

The goal of this paper is to use the ISS dimensional argument, for the first time in some examples, to derive the 
KZ spectra and the scaling of the interaction coefficients for a large number of commonly considered applications of 
weak turbulence. We then use our results to discuss the uniformity of the weak turbulence approximation in k for 
these physical systems. In the final section we consider the water wave system in more detail. It is shown that by 
considering the effect of both the gravity dominated and surface tension dominated parts of the spectrum together, 
one can build a consistent picture of energy transfer in the system, even when the flux is sufficiently large to cause 
breakdown of the weak turbulence approximation. 



II. DIMENSIONAL DERIVATION OF KOLMOGOROV-ZAKHAROV SPECTRA 



Before we begin, let us clarify a point of notation. We deal with isotropic systems. Physical quantities such as 
spectral wave- action density, nk, or spectral energy density, i?k = ^k^^k, only depend on the modulus, fc, of the 
wave-vector, k. It is often convenient to integrate over angles in k-space. We need to make a distinction between 
a spectral quantity which has been averaged over angle and one which has not. To do this, we use a regular type 
argument to denote a quantity which has been integrated over angles, as in n^, and a bold type argument to denote 
one which has not, as in nk. The two are easily related. Consider for example, the wave-action density: 

J Uk dk — y "-k dh 

= nDnkk°~\ (7) 

where fi/) is the solid angle coming from the integration over angles in _D-dimensional wave- vector space. We shall 
use C to denote a generic dimensionless constant whose value cannot be determined from dimensional arguments. 



A. Constant Energy Flux Spectra 

Suppose we have a wave system characterised by a single additional dimensional parameter. A, which appears in 
the dispersion relation in the form 

= Afc". (8) 

It is convenient to set the density of the medium to 1. Our unit of mass then has dimension L'^ and energy has 
dimension L^T~^. We suppose that the d-dimensional energy density, £, is finite in physical space. For example, 
d — 2 for water waves while d = 3 for acoustic waves. We denote the dimension of the Fourier transform used to go to 
a spectral description of the theory by D. Usually D = d but not always (see, for example, section rvT|) . The spectral 
energy density, Ek, is defined by 

Eudk= / Ekdk. (9) 
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Ek clearly has dimension L^~'^T~^. The energy flux, P, has dimension L^~''T~'^ and A has dimension L"T~-'^. Let us 
now consider the constant energy flux spectrum for this system. For 3-wave processes, the energy flux is proportional 
to the square of the spectral energy density so we can write 

Ek = CVPX^k^, (10) 

where C is a dimensionless constant and the exponents X and Y are to be determined by dimensional analysis. This 
yields 

X=i Y=i(d + a-7). (11) 

This argument, used by Kraichnan "i] in the context of Alfven waves, can be generalised to iV-wave systems. In a 
system dominated by N wave processes, the energy flux is proportional to the N — 1**^ power of the spectral energy 
density and a similar argument yields the scaling law 

Ek = CP^X^k^, (12) 

with X and Y given by 

2iV- 5 



X 



N - 1 



Y=(2a + d-6) + ^^. (13) 

Associated with each constant energy flux spectrum, we have a particle number (wave action) spectrum, rik. One can 
be obtained from the other via the relation 

) /-OO 

Ekdk = nD tJkn-kfc^^^rffc, (14) 
Jo 

where fl^) is the D-dimensional solid angle. The resulting scaling law for ?ik for an N wave system is 

nk = CP^A^/c^, (15) 

where 

X 



iV- 1' 

5 — Sa — d , , 

Y^ia + d-D-5)+ . (16) 



B. Constant Particle Flux Spectra 

In the case of a system with 4-wave interactions, the total particle number, N ~ J Uk dk is also a conserved 
quantity. As a result, there can also exist a constant-flux spectrum carrying a flux of particles rather than energy. 
Such behaviour is associated with a continuity equation of the form 

where Qk is the particle flux. One can perform the same dimensional analysis for this spectrum, bearing in mind that 
dimensionally, P = lui^Q. One obtains the following spectrum describing a constant flux of particles: 

or 

71k = CQUU-^(3D-2rf-a+10) (19) 
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C. Scaling of the Interaction Coefficients 

In the regime where the system is scale invariant, the nonhnear interaction coefficients, Vijk (3-wavc) and Tijki 
(4- wave) often possess nontrivial scahng properties. For the 3- wave case we have : 

V/(k/iki/(k2 = /I'^Vkkiks) (20) 

and for the 4-wave case : 

— ^''^Tkkik2k3- (21) 

In fact /3 and 7 cannot be arbitrary. They may be determined from dimensional analysis of the dynamical equations. 
Schematically, the dynamical equations for an N wave system look like 



dak 

dt 



+ Mk = j < '5(k) (dk)^-i. (22) 
Recalling that dimensionally, [k] — and [(5(k)] = L^, we see that: 



rr 1 ^'^^ 



< ttkO^/ >= S{k — k') rik- 



Determine the dimension of as follows, 



So 



M' = [(5(k-k')] 



6-dn 



This results in the following expression for the dimension of the interaction coefficient 

[i^] = T5(^-4)L5(^-2)(''-5), (23) 



and dimensional analysis then yields. 



where 



^iVki---k, 



A^(^-^)F«/k,...k., (24) 



7iv = -^{(iV-2)(d-5) + (Ar-4)a}. (25) 

Here /ki -kAr ^ dimcnsionless function of ki---kjv. Interestingly, for 4-wave systems, iV = 4, the scaling of 
the interaction coefficients depends only on the dimension, d, of the system and is independent of any dimensional 
parameter, including A. We see that all incompletely self-similar 4-wave systems exhibit the same scaling behaviour 
of their interaction coefficients, 

74 = 5 - d. (26) 
Applying our analysis to the 3- wave case yields iskikakg ~ k''^, where 

73 = i(5 + a-rf). (27) 

The scaling in this case depends on the dispersion index, a, but we see that 73 and a are not independent quantities. 
This fact, while obvious from this point of view, is possibly not fully appreciated. We shall see that the class of 
incompletely self-similar systems for which this analysis is valid includes most of the common physical applications of 
weak turbulence. 
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III. BREAKDOWN OF THE WEAK TURBULENCE APPROXIMATION 



By analysing the scaling behaviour of the kinetic equations describing nonlinear energy transfer in weak turbulence, 
Biven, Nazarenko and Newell have, under certain assumptions, given a set of criteria for the breakdown of the 
weak turbulence approximation. These assumptions are that the turbulent transfer is sufficiently local that after 
using homogeneity properties to remove the k-dependence of the collision integral and the other integrals arising in 
the expression for the frequency renormalisation, the remaining integrals converge. We discuss in the conclusion and 
in how this may not always be the case when the coefficient of long-wave short-wave interaction is too strong. 

In this section we shall adopt the commonly used notation 73 = /3 and 74 = 7. For three wave systems, breakdown 
occurs at small scales for /3 — 2a > and at large scales for /3 — 2a < 0. For four wave systems breakdown occurs at 
small scales for 7 — 3q; > in the presence of a pure energy flux and for 7 — 2a > in the case of a pure particle flux. 
The breakdown at large scales can be masked by the large scale forcing but breakdown at small scales, in the absence 
of dissipation, is taken to signal the onset of small scale intermittency. In the case of ISS systems we can construct a 
characteristic scale, /cnl, from the flux and the parameter A. From our previous discussion of dimensions, we see that 
the quantity, (^P/ X^y/^^'^^^'^ has the dimension of a length. Thus in the case of a finite energy flux, we define 



fcNL - ( ^ ) . (28) 



Likewise, in the case of a finite particle number flux, Q, in a four wave system we define 



ANL=(^j (29) 

For small fluxes, P ^ 0, we see that the breakdown occurs at small scales for 5 — d — 3a > for finite P and for 
5 — d — 2a > for finite Q. Upon substitution of (|23(l into these expressions, we recover the criteria of 2] in terms of 
the scaling exponents of the interaction coefficients. It is interesting to note that for finite energy flux, the breakdown 
criterion is a < 2/3 in 3 dimensions and a < 1 in 2 dimensions. However it is known from the work of Krasitskii 
that for a < 1, 3- wave terms in the interaction Hamiltonian are nonresonant and can be removed by an appropriate 
change of canonical variables to give an effective description in terms of four wave interactions. Thus the small scale 
breakdown criterion can never be realised for three wave systems in two or three dimensions. This means that a 
significant number of physical systems cannot be hoped to exhibit intermittency at small scales as discussed in the 
following section. Conversely, these systems can always exhibit intermittency at large scales, provided the forcing is 
sufficiently strong, without affecting the validity of the weak turbulence approximation at small scales. At this point, 
it is worth mentioning that the case a = 1 is borderline in two dimensions. Such three wave systems, 2-d sound being 
an example, are known to be rather special and must be carefully treated separately. 



IV. EXAMPLES OF 3- WAVE SYSTEMS 

Sound and magnetic sound in 3 dimensions 

Acoustic turbulence 0, 0, corresponds to the almost linear dispersion Wk ~ ck, where c is the sound speed or 
magnetic sound speed, so that a = 1 and d = 3. We thus obtain the following pair of spectra for the energy and wave 
action 

Ek = C^/Vck-i , nk = C'\-k-^. (30) 

V c 

These are the original spectra obtained by Zakharov and Sagdeev. According to our analysis, this spectrum remains 
uniformly valid at small scales. 



3-D Alfven waves 



3-D Alfven wave turbulence was originally considered by Iroshnikov |lO| and Kraichnan (J] in the 60 's . Such waves 
are also weakly dispersive and from the point of view of the dimensional analysis are identical to acoustic waves 
discussed above. The resulting — | and — | spectra are not actually realised in real plasmas because the true Alfven 
wave turbulence is anisotropic. 
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Quasi 2-D Alfven waves 



In reality, the Alfven turbulence is strongly anisotropic and is described by quasi-2D rather than 3D spectra. For 
this system we have again a — I but d = 2. This yields the stationary spectra 

Ek = Cy/P^k-^, rik = C'\-k-^, (31) 

V c 

which are the spectra obtained using a dimensional analysis by Ng and Bhatachargee and analytically derived by 
Galtier, Nazarenko and Newell As mentioned already, the case a = 1 is borderline in 2-d so that our argument 
concerning the breakdown of the weak turbulence approximation remains inconclusive. In fact, unlike typical three 
wave systems, this system does exhibit breakdown at small scales as shown in |l2| . 

Capillary waves on deep water 

A — where a is the coefficient of surface tension. In this case, a — ^, d — 2 and the Kolmogorov spectrum is 

Ek^C^/Pah-^, ni, = C"^/Pa-ik-^. (32) 

This spectrum was first derived by Zakharov and Filonenko . There is no small scale intermittency in this system. 



V. EXAMPLES OF 4- WAVE SYSTEMS 



Gravity waves on deep water 

For this system, a ~ ^, d = 2 and A = where g is the gravitational constant. The Kolmogorov spectrum 
corresponding to a constant flux of energy is then 

Ek^CP^gh-i, nu^C'Ph-^. (33) 

There is also a second spectrum corresponding to a constant flux of wave action, 

n^^CQ^gh~'f. (34) 

These spectra were obtained by Zakharov and Filonenko |l4j . In this case, the energy spectrum breaks down at small 
scales. 

This is one of the cases where certain integrals appearing in the frequency renormalisation series diverge on the 
K-Z spectrum. The problem is that the interaction coefficient between the high k mode and long wave partners in its 
resonant quartet is too strong. This leads to a modification of the breakdown criterion and means that the breakdown 
can occur for values of k less than that value calculated when local interactions dominate. 



Langmuir waves in isotropic plasmas, spin waves 



Langmuir waves are described by the dispersion relation 

u;l - (1 + 3rlk^) , (35) 

where ujp and rjj are the plasma frequency and Debye length respectively. Magnetic spin waves in solids also obey a 
dispersion relation of this type but the physical meaning of the dimensional parameters is different. For long 
Langmuir waves we can Taylor expand Wk as uj-k, = ujp + ^ajprj^k'^. The constant factor cUp cancels out of both sides 
of the 4- wave resonance condition so that the effective dispersion is ~ fc^. Thus taking A = ^pf'o^ « = 2, d = 3 we 
obtain the energy spectrum 

Ek^CP^ {LUprl,)k-^, nk = C' pi k-^. (36) 
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FIG. 1: Schematic representation of the different scaling regimes in the surface wave problem. The width of the window of 
fully nonlinear turbulence depends on the energy flux input at the forcing scale. 



Using equation it^ . /cnl {P / uj^r^^)^ / ^ — > as P ^ so this spectrum should break down at large scales. The 
second spectrum carrying the wave action flux is 



which should also break down at large scales since /cnl ~ (Q/'^p'"i>)^^^ 
spectra were originally derived by Zakharov 1^ (1972). 



(37) 

as Q ^ from equation (|25J). These 



VI. A 5- WAVE EXAMPLE 



One-dimensional gravity waves were considered by Dyachenko et al 17]. They found that the 4- wave interaction 
coefficient is identically zero on the resonant manifolds so that the nonlinear exchange of energy in the system is, in 
fact, due to 5- wave interactions. This system is an example of a case where the dimension of the physical energy 
density differs from the dimension of the Fourier space. In this case, d = 2 but D — 1. For this system, A = and 
a = i. Applying H15|) and (|12|l we obtain the finite energy flux spectrum 

Ek = CPigh^, (38) 
nk = C'{P^)-^k^, 

as found in From equation (|28|l . /cnl ^ P '^^^g goes to large k for small P. Therefore this spectrum should 

break down at small scales. 



VII. MATCHING THE GRAVITY AND CAPILLARY WAVE SPECTRA 



Let us now consider the complete surface wave problem including both gravity and surface tension effects. At large 
scales the system is entirely gravity dominated. We assume that the forcing is at large scales only. At small scales 
the system is entirely surface tension dominated down to the viscous scale where the wave energy is finally dissipated. 
The characteristic scale, fco, where surface tension and gravity are comparable can be estimated from the dispersion 
relation, 

uj{k) ^ ^/gk + ak^. (39) 
The gravity and surface tension effects are of comparable order when 



(40) 
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We expect that for k << ko the system is well described by the gravity wave spectrum (|33|) . for k >> kg the system is 
well described by the capillary wave spectrum, 132|l . In between there is a non-scale invariant cross-over regime. Let 
us consider the question of whether the weak turbulence approximation remains consistent through this cross-over 
regime. In order for the turbulence to remain weak as we approach kg from the left, the gravity wave spectrum must 
remain valid at least to the scale ko where surface tension effects can start carrying the flux. Thus we require 



4l > fco (41) 



where /c^l breakdown scale for pure gravity waves which we calculate from (|28|l : 

« P-h- (42) 
Using expressions (|42|) and H4U|I this gives us a condition on the flux, 

P<iga)-^ (43) 

In order for the turbulence to remain weak as we approach fco from the right, the capillary wave spectrum should be 
valid by the time we reach scale ko so that it can connect to the gravity wave spectrum. Thus we require 

4t < ko (44) 

where k'^^ is the breakdown scale for pure capillary waves, 



kl,i^pia-\ (45) 



Inserting expressions and (gnj) this gives us the same condition, l|43|) . on the flux! We see that there is a critical 
energy flux, Pc = {gfj f^^ which can be carried by the weak turbulence spectra. The issue of what happens if P > Pc 

is of paramount interest. It is clear that in this case there is a window in k space corresponding roughly to fc^L' ^nl 
where the nonlinearity is not weak and the dynamics is presumably dominated by fully nonlinear structures. This 
situation is illustrated schematically in figure^ It is suggestive that the value of Pc, if expressed in terms of the wind 
speed, corresponds roughly to the threshold for the formation of whitecaps on the ocean surface [isj. 

The phenomenon of intermittency is thought to be associated with the generation of such strongly nonlinear 
structures and would manifest itself in a deviation of the structure functions, S'Ar(ri, . . . r^r-i), of the wave field from 
joint Gaussianity. If one assumes that the statistics are dominated by whitecaps then one can estimate the scaling 
behaviour of field gradients. However, the support of the set of singularities need not be simple set. There are reasons 
to expect that whitecaps are supported on a fractal set of dimension < D < lHJj although we will consider fractal 
sets up to dimension 2. In this case the rt*'* moment of the field gradients scales as 



Sj^ir) ^ {Aer [^) , (46) 

where A9 is a characteristic size of the jump discontinuities in the derivative and L is the integral scale. It then 
follows (see 

13 sec. 8.5) that, 

S2Nir) /r\(i-^')(2-o) 



(^2(r)) 



N 



(z) • (4^) 



For D <2 the system deviates from joint Gaussianity. Such behaviour is generally thought to be beyond the standard 
picture of weak turbulence. Nonetheless, Biven, Nazarenko and Newell |3, 1^ have calculated the first correction to 
joint Gaussian statistics in the case where the weak turbulence approximation breaks down. For TV even, 



S2N{r) 
(52(r))^ 



N/2 



^C^. (pi/V"-i) " +... (48) 



Breakdown occurs for 7 > 3a in which case, the second term in (|48|l scales like r^°'~''/^)'^^~^) as r ^ 0. If we wish to 
attribute this breakdown to the emergence of whitecap-dominated behaviour, we observe that it is possible to match 
the scalings ()47|l and H48() for all N if we choose 

D = a-^ + 2. (49) 
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For gravity waves, a = 1/2 and 7 = 3 so the dimension of the set of whitecaps would be 3/2. It would be nice if 
one obtained a value of D less than one but there are several reasons why this argument is an oversimplification. In 
particular, expression 148|l represents only the first terms in an infinite series. It is highly likely that in the regime 
where weak turbulence breaks down, the higher order corrections which are neglected here actually contribute strongly. 
Nonetheless it is a nontrivial fact that this matching can be done consistently for all values of N simultaneously, even 
if the actual value of the fractal dimension obtained here must be considered with caution. 



VIII. SUMMARY AND CONCLUSION 



Our aim in this article was to show that many of the commonly considered applications of weak turbulence possess 
the incomplete self-similarity property, which can be exploited to obtain core results using a simple dimensional 
argument without resorting to the more complex methods required in general. For such systems, recent results on 
the breakdown and range of applicability of weak turbulence can also be obtained in a simple way. It was found 
that dimensional considerations rule out the development of small scale intermittency in most physically relevant 
three- wave systems. 

We considered the gravity-capillary surface wave system in more detail and discussed, from the point of view of the 
dimensional quantities present, how the validity of the weak turbulence approximation depends only on the energy 
flux input at the largest scales. Even in the case where this flux is large enough to cause breakdown of the gravity 
dominated part of the spectrum, we still have a consistent mechanism to transfer energy to the viscous scale which 
consists of two different weakly turbulent regimes connected by a window of scales where the nonlinearity is strong. 
We made some speculative observations about the relationship between the breakdown of weak turbulence and the 
emergence of whitecaps in this window of strong turbulence. 

It is appropriate that we finish with some balancing remarks about situations in which the simple approach outlined 
here does not work. Firstly, there are obviously cases of physical interest which are self-similar. Some important 
examples are provided by optical waves of diffraction in nonlinear dielectrics and the turbulence of waves on Bose- 
Einstein condensates, both of which are described by the Nonlinear Schrodinger equation j2^. In these cases, there 
are two relevant dimensional parameters. 

Secondly, even in the case of incompletely self-similar systems, a cautionary note should be sounded. The long time 
behaviour of these systems is determined by the kinetic equation for the spectral wave action density, 

^=T2K]+T4K]-f (50) 
at 

and a nonlinear frequency modulation, 

Wfe — > Wfc -I- VI2 [uk] + ^4. [nk] + (51) 

For 4- wave systems with an interaction coefficient, Tkkik2ka, 

T2[nk] = 0, (52) 
Ti [uk] = J \Tkk1k2k3 f "i'T-2«3<5(^^oi,23)<5(koi,23)dki23 -l- rifelm fli [uk] , (53) 

with the frequency modulation integrals given by 



^2[nk] = J Tkkikki-rii dki, (54) 
lmV,4[nk] = / iTfefc^fe^fegl^ (711713 4- 711712 " n2n3)(5(tJoi,23)<5(koi,23)rfki23. (55) 



Notice that the imaginary part of the frequency modulation enters the collision integral. It has been pointed out 
by Zakharov and others that if the interaction coefficient is not uniformly homogeneous in its arguments, 

Tkk,k2k, ~ {kk2)^^-^'>''{hk3)^^ < X < 1, (56) 

for example, the frequency correction, Im ^2 [nk] can be divergent at low k, even though T is still homogeneous of 
the same degree, 7. Luckily, this divergence cancels in the kinetic equation to lowest order, T4 [uk]- However, the 
divergences in the fi's may resurface at higher orders in the full collision integral, H5U|) . It is not clear whether we are 
rescued by such cancellations as occur in T4 [nk]. 
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If these divergences persist then the dimensional argument appHed in section lllll to estimate the range of vaHdity 
and breakdown scales would be not valid and would require modification to include some non-universal dependence 
on X and a low k cutoff. Thus wave turbulence may possess a mechanism for breaking incomplete self-similarity even 
without the need to introduce additional dimensional parameters. Whether this occurs in reality is, at present, an 
open question. 
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